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A simulation of the Neolithic transition in the Indus valley 

Carsten Lemmen and Aurangzeb Khan 



The Indus Valley Civilization (IVC) was one of the first 
great civilizations in prehistory. This bronze age civilization 
flourished from the end of the fourth millennium BC. It dis- 
integrated during the second millennium BC; despite much 
research effort, this decline is not well understood. Less re- 
search has been devoted to the emergence of the IVC, which 
shows continuous cultural precursors since at least the sev- 
enth millennium BC. To understand the decline, wc believe 
it is necessary to investigate the rise of the IVC, i.e., the es- 
tablishment of agriculture and livestock, dense populations 
and technological developments 7000 3000 BC. Although 
much archaeological information is available, our capabil- 
ity to investigate the system is hindered by poorly resolved 
chronology, and by a lack of field work in the intermedi- 
ate areas between the Indus valley and Mesopotamia. We 
thus employ a complementary numerical simulation to de- 
velop a consistent picture of technology, agropastoralism and 
population developments in the IVC domain. Results from 
this Global Land Use and technological Evolution Simulator 
show that there is (1) fair agreement between the simulated 
timing of the agricultural transition and radiocarbon dates 
from early agricultural sites, but the transition is simulated 
first in India then Pakistan; (2) an independent agropas- 
toralism developing on the Indian subcontinent; and (3) a 
positive relationship between archeological artifact richness 
and simulated population density which remains to be quan- 
tified. 



1. Introduction 

The Indus Valley Civilization (IVC) , often termed Harap- 
pam civilization after its major type site, flourished along the 
banks of the river Indus and its tributaries, including the ad- 
jacent coastal areas, between the fourth and second millen- 
nium BC. The IVC is characterized by urban centers, bronze 
technology and seals, trade networks with Mesopotamia and 
Arabia, and an as yet undcciphcred writing system [Shaffer, 
1992; Possehl, 1998; Kenoyer, 2008]. Much research has 
been devoted to the disintegration of the IVC during the 
second millennium BC, and much less to the emergence of 
this great civilization which shows continuous cultural pre- 
cursors at least since the seventh millennium BC, such that 
some authors prefer to speak of the Indus valley cultural 
tradition [e.g., Kenoyer, 2006]. 

The Indus valley is one of the two great river basins on 
the Indian subcontinent, separated from the Ganges val- 
ley to the cast by the Aravalli mountain range and the 
Thar desert (Figure lb). The region has been occupied by 
anatomically modern humans for at least 34,000, but pos- 
sibly >80,000 years [Petraglia et al, 2010]. The late Pa- 
leolithic and the Mesolithic are visible in stone and blade 
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industries; microlith blades indicate composite tools, differ- 
entiated occupation levels and a broad-spectrum diet of for- 
aging people from 26,000 years before present (a BP) until 
at least the second millennium BC, well into the Bronze age 

[Petraglia et al, 2009]. 

Neolithic subsistence is based on a combination of bar- 
ley and wheat dominated agriculture and herding of sheep, 
goats, and cattle. Domesticates are of mixed origin: barley 
and wheat have no local wild progenitor, but zebu cattle 
bears a clear signature of local domestication [Chen et al, 
2010; Patler, 2011a]. The beginning of the Neolithic on the 
Indian subcontinent is regionally diverse and ranges from 
6500 BC in Baluchistan [Jarrige, 1995] to 3000 BC or later 
on the central Indian plateau and in the Himalaya mountain 
range [Samuel, 2000; Boivin et al, 2008]. Our focus here is 
on the period which spans from the Neolithic to the begin- 
ning of the Bronze age of the IVC (7000-3200 BC), i.e., on 
the early food production and regionalization eras, or the 
Mehrgarh phase and parts of the early Harappan phase ac- 
cording to the different chronologies [Sharma et al., 1980; 
Fuller, 2006; Kenoyer, 2008]. 

Possehl [2002] separated these phases and eras of the 
pre-Harappan Neolithic into two stages (Figure la): Early 
Neolithic and Developed Neolithic. The Early Neolithic 
is divided in the Kill Ghul Mohammad (KGM) and the 
Burj Basket-Marked (BBM) periods. Sites attributed to 
these periods are found in the western borderlands of the 
greater Indus region and cluster in the hills and piedmont of 
Baluchistan (Figure 1 bl b2). During the BBM, soft pottery 
and handmade pottery are introduced. Deposits of trash 
with burnt pebbles, ash, animal bones, bone tools, hammer 
stones, polishers colored with red ochre, and a large collec- 
tion of blades, cores, and flint debris have been discovered, 
and point to leather making, basket making, or weaving. 
Domesticated wheat, barley and cotton were used, pulses 
may have been present; the use of copper is visible in the 
form of arrowheads and beads [Moulherat et al, 2002; Fuller, 
2011a]. 

The Developed Neolithic is characterized by improved 
pottery, developed agropastoral communities, growth, con- 
tinuity, and geographical expansion. During its first period 
(here defined as both Togau and Sheri Khan Tarakai, SKT, 
Figure 1 b3), the settlement patterns show growth of village 
life; technological innovations include the use of gold, the 
manufacture of compartmented seals, glazed steatite, and 
beads. An eastward expansion is visible in sites emerg- 
ing along the ancient Ghaggar Hakra river. The second 
period of the Developed Neolithic (Figure 1 b4) comprises 
several cultural complexes. In the Kechi Beg complex in 
Baluchistan wheat replaces barley; an irrigation canal near 
Mehrgarh may date to this period [Kenoyer, 2008]. On the 
Punjab plain, a change in the assemblage of pottery types 
characterizes the Hakra Ware complex: Microlithic tools are 
abundant, and sites represented by a light scatter of pottery 
without a buildup of a midden point to seasonal occupations 
[Possehl, 2002]. At the same time, the Anarta complex ap- 
pears in Gujarat with distinct pottery styles. 

Our capability to investigate the system more closely and 
consistently in the regional context is hindered by the poorly 
resolved chronology, and by a lack of field work in the inter- 
mediate areas between the Indus valley and Mesopotamia. 
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Figure 1. Temporal and spatial domain for this study: a) Neolithic chronology of Baluchistan, the Indus valley, and 
Gujarat based on Possehl [2002] with updates from Fuller [2006]. Relative resource shortage is based on palaeoproxy 
evidence from Lemmen and Wirtz [2012] and a transient Holocene climate simulation by Brovkin [2002]. b-e) Geography, 
topography and site distribution of artifacts [Law, 2008] typologically associated with the Neolithic. 



We thus employ a complementary numerical simulation to 
the IVC domain. Ackland et al. [2007] simulated the Ne- 
olithic transition in the IVC with an advancing frontier 
spreading from Southwest Asia, and a competition between 
migrating original farmers and resident converts. After 5000 
simulation years (9000-4000 BC), their simulated spread of 
agriculture into the Indian subcontinent exhibits a demar- 
cation line along the Indus river, separating original farmers 
to the west from converts to the east. 

A local South Asian agricultural center was considered 
in the numerical model by Patterson et al. [2010]. They 
assumed a background hunting-gathering population, and 
initial farmer populations near Mehrgarh and at the eastern 
edge of their model domain. After 2000 simulation years, 
the population in the Indus valley is dominated by expand- 
ing original farmers, a lesser share of converts, and a tiny 
fraction of remaining hunter-gatherers. 

We go beyond the biophysical approaches by Ackland 
et al. [2007] or Patterson et al. [2010]; in our Global Land 
Use and Technological Evolution Simulator [GLUES, e.g., 
Lemmen et al., 2011], we include socio-technological innova- 
tion in addition to migration, population growth and subsis- 
tence change to develop a consistent picture of the dynamics 
of technology, agropastoralism and population developments 
during the Mesolithic Neolithic transition. This simulated 
data is subsequently compared to a representative archae- 
ological data set [Law, 2008], where the chronology of Ne- 
olithic sites is used to detect the regional transition from for- 
aging to farming, and where the site occurrence frequency 
is used to detect the Neolithic demographic transition. 



2. Data and model description 

2.1. Global Land Use and technological Evolution 
Simulator 

GLUES is a socio-technological model which hindcasts 
technological evolution, potential population density, and 
the timing of the transition to agropastoralism for 685 world 
regions based on the geoenvironmental and cultural con- 
texts, and innovation and adaptation of regional popula- 
tions. Each local population in an approximately country- 
size region utilizes its natural resources and interacts with 
its geographical neighbors via trade and migration. The full 
model is described by Wirtz and Lemmen [2003] , with refine- 
ments of geography [Lemmen, 2009], migration and knowl- 
edge loss [Lemmen and Wirtz, 2010], and climate events 
[Lemmen and Wirtz, 2012]; it has boon validated against 
archaeological data by Lemmen et al. [2011]. Key concepts 
arc summarized below: 

Local societies are defined by their population density 
(P) and by three characteristic traits: (1) technology (T), 
which describes efficiency of food procurement; (2) share 
of agropastoral activities (Q) as the allocation pattern of 
manpower to farming and herding as opposed to foraging; 
and (3) economic diversity (A*"), a technology-related trait 
which represents the different subsistence styles of a com- 
munity. The temporal change of traits follows the direction 
of increased benefit for success (i.e. growth) of its associ- 
ated population with a speed related to the trait's variance 
[Fisher, 1930; Kisdi and Geritz, 2010]. Population growth 
rate (r) depends on resource availability (E), subsistence 
intensity, overexploitation of resources, administrative over- 
head and health improvements; the latter four terms are 
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functions of P, T, N, and Q, while the resource availabil- 
ity E is provided externally. E is estimated from the net 
primary productivity, which is in turn based on Climber-2 
[Claussen et al., 1999] transient temperature and precipi- 
tation anomalies from the IIASA climatological data base 
[International Institute for Applied Systems Analysis, Lee- 
mans and Cramer, 1991]. 

Abrupt climate deteriorations modulate E; they are in- 
corporated by using the extreme event anomalies found 
in 124 globally distributed palaeoclimate proxies collected 
and analyzed by Wirtz et al. [2010]; the estimation of re- 
source failure from proxy climate is described by Lemmen 
and Wirtz [2012]. Figure la shows the average relative de- 
cline of E in the Indus valley region; these anomalies are 
derived from regional abrupt changes found in oxygen iso- 
tope anomalies from Kyrgyzstan [Ricketts et al, 2001], Ti- 
bet [Pontes et al, 1996], the Arabian Sea [Doose-Rolinski 
et al., 2001; Gupta et al., 2003]. and Oman [Fleitmann et al., 
2007]. 

We set up the eight global model parameters and initial 
values such that the simulation is able to hindcast an ac- 
curate timing and location of the early farming centers in 
Southwest Asia, northern China, and Mesoamerica [Smith, 
1998], and a reasonable global pattern of the subsequent Ne- 
olithic transition. The simulation is started at 9500simBC 
(simulation years BC; equivalent to calendar years EC, 
but emphasizes the artificial model time scale, Lemmen 
et al. 2011). All of the 685 biogeographically defined re- 
gions, including 32 around the Indus valley, are initially 
set with farming activity at 4%, 1/4 established agropas- 
toral economies and unit technology; this setup represents a 
low density Mesolithic population relying on a broad spec- 
trum foraging lifestyle with low unintentional farming activ- 
ity. The regional simulated transition to agropastoralism is 
recorded when Q exceeds 0.5. 

2.2. Artifact database 

We used the Indus Google Earth Gazetteer (version Au- 
gust 2008) compiled by Law [2008] for the geolocation of 
artifacts relating to the Indus valley. From this database 
of 21687 metal and stone artifacts, we use here the cultural 
attribution and location of those 2028 artifacts and 368 sites 
that are in the spatial and temporal domain of our study. 
For KGM, 18 sites contribute 90 artifacts to the data base; 
for BBM, there are 159 artifacts at 32 sites. We combined 
the sites attributed to the Togau period with sites from the 
SKT (both with starting date 4300 BC, Fuller 2006) and 
obtain 317 artifacts from 65 sites. 

The subsequent period is represented by the Hakra Ware, 
Kechi Beg, and Anarta complexes. Combined, there are 
1250 artifacts (396, 670, and 184, respectively) from 253 sites 
(99, 112, and 42, respectively). While the main chronology 
is based on Possehl [2002], sites from the Anarta complex 
were included in this period based on Fuller [2006]. To our 
knowledge, this gazetteer is the most representative data 
set of lithic and metal artifacts of the Indus valley tradition. 
It has been used to investigate the trade and distribution 
networks of the IVC by its author [Law, 2011]. 

In archaeological or palaeoenvironmental archives, the de- 
crease of the number of dates further back in time that 
is caused by gradual destruction of the archive by post- 
depositional processes results in a taphonomic bias [Surov- 
ell et al., 2009]. We applied to the artifact density data 
a bias correction derived from the comparison of volcanic 
eruption layers in the GISP2 (Greenland, no taphonomic 
bias) ice core with radiocarbon dated ash deposits; the 
correction factor c depends on time t (in a BP): c{t) = 
5.73-10^(2175.4-|-t)-^ * [Surovell et al., 2009]. This method 
was successfully validated with Canadian coin mintage data 
compared to a coin collection and applied to North Ameri- 
can radiocarbon dates by Peros et al. [2010]. 



3. Results and Discussion 

The simulated transition to predominantly agropastoralism- 
based subsistence occurs in the spatial domain of the 
Indus valley and surrounding areas between 6300 and 
3800 simBC (Figure 2a). There is a general east- west trend 
in the simulated transition dates. Contemporary with the 
KGM, the earliest transitions are simulated in India (bo- 
fore 5600 simBC), including Rajasthan, Gujarat, and the 
Ganges valley, and early transitions (before 5000 simBC) 
in southern Sindh and northern Pakistan, including parts of 
the Punjab and Kashmir. During this first period, 15 out 
of 32 simulation regions undergo the transition to agropas- 
toralism. Most of the remaining regions (Pakistan, eastern 
Afghanistan, eastern Iran and Makran coast) transition in 
the BBM period before 4300 simBC. Of these, agropastoral- 
ism is simulated first in the Punjab and northern Baluchis- 
tan (before 4700 simBC), then in northern Sindh (before 
4500 simBC). The simulated transition occurs latest, dur- 
ing the Togau period, in the eastern model domain, from 
central Afghanistan to the Pakistan-Afghanistan-Iran bor- 
der triangle (4300-3800 simBC). 

A detailed and quantitative comparison is limited fore- 
most by the inherent timing uncertainties in the typolog- 
ically dated artifacts. This is closely related to the sec- 
ondary comparison difficulty, which is created by the lack of 
temporal continuity of samples from the Mesolithic to the 
Neolithic. While we have shown how these issues can be ad- 
dressed in the context of plenty and precise dates for Europe 
[Lemmen et al., 2011], the comparison with dates from the 
Indian subcontinent can only be qualitative at this moment 
(see also Patterson et al. 2010, who make a qualitative com- 
parison not with single sites, but with approximate timing 
contours). 

In addition, the coexistence of Mesolithic and Neolithic 
subsistence in the archaeological record complicates the 
comparison between model and data, because the degree of 
both subsistence styles carmot be quantified with the current 
data. This coexistence is evident in the long-term contin- 
ued use of (Mesolithic) microlithic tools, transhumance life 
style, i.e., the seasonal back-and- forth movement of herders 
with their livestock, and the continuation of foraging prac- 
tice weU into the Neolithic [Mughal, 1990; Fuller, 2011a], it 
is not resolved in the model, where agropastoral life style is 
usually dominant within less than 500 years after its emer- 
gence [Lemmen and Wirtz, 2012]. 

The model simulates first agropastoralism in northern 
Baluchistan and southern Sindh in the KGM period consis- 
tent with earliest sites in these areas. For the major part of 
Baluchistan, including the Mehrgarh site, the model places 
the transition in the BBM period, which is later than the 
earliest site dates, but agrees with dates from many other 
locations in Baluchistan. The most obvious mismatch is the 
early simulated (KGM and BBM) transition along the In- 
dus, the Ghaggar-Hakra and Punjab rivers, where site dates 
are absent until the Kechi Beg and Hakra periods. 

While evidence for early archaeological sites may have 
been lost in the geomorphologically dynamic Indus river 
floodplain, the lack of early sites in Rajasthan is striking, 
considering that the model sees a biogeographically favor- 
able environment and cultural setting for an early Neolithic: 
one should ask whether an agricultural center on the Indo- 
Gangetic divide is undetected so far. First evidence for 
northern Indian agriculture could come from further east. 
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Figure 2. a) Timing of the Mesolithic-Neolithic transition from a GLUES simulation (background color, simulation 
years BC) and inferred from the presence of Neolithic sites (foreground triangles, calendar years BC) in the greater Indus 
valley, b) Recovered artifact number in relation to simulated population size. Histogram of artifact density in four pre- 
Harappan periods of the Indus valley and Baluchistan based on Law's (2008) Indus Google Earth Gazetteer, and corrected 
for taphonomic bias (bars, annotated with n x c, where n is the raw number of artifacts and c the taphonomic correction 
factor); GLUES-simulated population size in all areas within 200km from where Neolithic pre-Harappan artifacts were 
discovered (solid line). 



from Lahuredawa in the middle Ganges plains, where ev- 
idence for domesticated rice appears as early as 6500 BC 
[Tewari et al, 2008] (but see Fuller 2011a for a debate on 
the rice morphology) . Less favored by the model are the val- 
leys along the Indo-Iranian plateau, where broad subsistence 
possibilities are seen as one precondition for the rise of the 
IVC and where agropastoralism arose before 6500 BC [Jar- 
rige, 1995]. The model might underestimate the potential 
for agropastoralism in this area because of its coarse spatial 
scale. 

From a regional perspective, the simulated east-west di- 
rection of agropastoralism contradicts the eastward trend 
seen in the site data. The westward direction had been 
suggested in earlier literature on the Indus valley tradition: 
Wheeler [1959] wrote of "movement of the Neolithic from 
Burma if not behind." In their simulation of the spread 
of agriculture into the Indian continent, Patterson et al. 
[2010] assume such a Southeast Asian center in addition 
the Mehrgarh location. While Wheeler's (1959) earliest 
dates (1000 BC) have been pushed further back by more 
recent datings, this westward view cannot be upheld from 
either artifacts [Law, 2008] or the description of sites in cen- 
tral India (3000 BC, Boivin et al. 2008), the Indus plain 
(4000 BC, Mughal 1990), northern India (5000 BC, Sharma 
et al. 1980), or Baluchistan (6500 BC, Jarnge 1995); all of 
these demonstrate the opposite east- and then southward 
trend. 

The simulated earliest agropastoralism in the Indus val- 
ley region is independent of the two prescribed early cen- 
ters of agropastoralism in Eurasia, i.e., China and South- 
west Asia; towards the end of KGM, however, a continu- 
ous band of early agropastoralism ranges from the north- 
ern Indian subcontinent through northern Southeast Asia 
to China (not shown). This South and East Asian agricul- 
tural area connects to the Southwest Asian center during 
the BBM period. The simulated independent emergence of 
South Asian agropastoralism corresponds well to evidence 
for independent domestication of cattle and rice [Chen et al, 



2010; Fuller, 2011b]. Though not completely detached from 
Southwest Asian domesticates assemblage in the data, this 
separation is visible in the delayed simulated transition to 
agropastoralism in the area between the Indo-Gangetic and 
the Southwest Asian founder centers. 

Building on the earlier idea of Sanders et al. [1979], we 
relate the number of artifacts in the database (corrected for 
taphonomic bias) to population estimates from our simu- 
lation (Figure 2b). We summed the total simulated pop- 
ulation in each half-degree model grid cell within 200 km 
distance to the find sites and obtain a population trajectory 
which shows an increase from 0.2 million at 6000 simBC to 
2.5 million at 3500 simBC; this corresponds to a population 
density of 0.1 and 1.3km^^, respectively. At 3500 simBC 
our estimate is higher than the one by Patterson et al. [2010, 
0.2 km~^], as should be expected from our calculation from 
the surroundings (within 200 km) of mostly settlement sites. 
To first approximation, the increases in simulated popula- 
tion size and in reconstructed artifact density correspond: 
the model correctly represents the emergence of dense set- 
tlements. Statistically, however, no meaningful relationship 
should be derived from only four different periods. 

In summary, there is fair qualitative agreement be- 
tween the reconstructed and the simulated transition to 
agropastoralism, both in the spatio-temporal distribution of 
agropastoral sites and simulation regions, as well as the over- 
all artifact and population density. Important lessons from 
the model are that it corroborates an independent agropas- 
toral center, and that it suggests potential for earlier Indian 
agropastoralism. It will be left to further studies to quan- 
tify the agreement and to elaborate on the reasons, both in 
archaeology and in the simulation, of discrepancies. 

4. Conclusion 

We presented a novel numerical simulation of the pre- 
Harappan Neolithic of the Indus valley tradition in the con- 
text of representative radiocarbon dates from material cul- 
ture. Simulated population size is qualitatively reflected in 
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artifact frequency in four pre-Harappan periods. Within the 
uncertainties of the coarse chronology, the model predicts 
the spatio-temporal structure of the Neolithic transition of 
this area fairly well. Our simulation points to a possible 
earlier Neolithic in northern India than seen in the data, 
and it corroborates an independent South Asian Neolithic. 
For quantitative model-data comparisons and clarification 
of the connection to Southwest Asia, better chronologic con- 
trol of the pre-Harappan material is needed as well as more 
evidence from the intervening Iranian plateau. 

Supplementary material 

The simulated data have been permanently archived on 
and are freely accessible from PANGAEA (Data Publisher 

for Earth & Eiivirormiental Science, http://www.pangaoa.de) 
as a notCDF datasct with reference "Lemmen, C. and Khan, 
A. (2012): Simulated transition to agropastoralism in the In- 
dus valley 7500-3000 BC. Dataset #779706, doi:XXXX". 
GLUES is open source software and can be downloaded from 
http:/ /glues. sourceforge.net. 
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